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AFIT/GAK/AA/81D-18 


Abstract 


Modern  optimal  control  techniques  are  used  to  design 
a  digital  controller  for  a  large,  flexible  space  structure. 
The  structure  is  modeled  as  a  tetrahedron  formed  by  four 
lumped  masses  connected  by  massless  truss  members.  A 
NASTRAN  analysis  is  used  to  generate  the  twelve  mode  shapes 
and  frequencies  of  oscillation. 

The  controller  is  designed  in  discrete  time  using 
linear  optimal  regulator  theory.  The  states  consist  of 
mode  amplitudes  and  velocities.  Because  these  amplitudes 
and  velocities  are  not  directly  available  to  the  controller, 
an  observer  is  used  to  estimate  the  states.  The  feedback 
gains  for  both  controller  and  observer  are  determed  using 
steady-state  optimal  regulator  theory.  The  controller  is 
designed  based  upon  only  four  modes  of  the  system.  The 
four  highest  frequency  modes  are  truncated  from  the  design 
to  signify  system  reduction. 

Active  control  is  achieved  using  six  collocated 
sensor-actuator  pairs.  The  effectiveness  of  the  controller 
is  determined  using  the  pointing  accuracy  as  a  figure  of 
merit.  The  controller  is  judged  by  its  ability  to  bring 
the  "pointer"  back  .o  the  nominal  line  of  sight  from  an 
initial  displacement. 

Three  controller  designs  are  compared.  The  first 
demonstrates  the  effects  of  "control"  and  "observation 
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spillover".  The  second  uses  ;i  singular  value  decomposition 
technique  to  generate  a  gain  matrix  to  eliminate  the  spill¬ 
over  terms.  The  second  design  also  uses  a  simplified  dis¬ 
crete  performance  index  for  optimal  control.  The  third 
design  uses  the  same  technique  as  the  second  to  eliminate 
the  control  and  observer  spillover,  but  it  uses  a  more 
accurate  performance  index. 

The  first  design  proved  to  be  inadequate  to  achieve 
the  desired  performance  because  of  the  spillover  terms.  When 
the  gain  was  increased  to  achieve  the  desired  performance, 
the  system  became  unstable.  Both  designs  II  and  III  were 
able  to  achieve  the  desired  performance  for  small  sampling 
times.  For  increasingly  larger  sampling  times,  system  per¬ 
formance  was  degraded  though  for  both  designs  the  system 
remained  stable.  In  addition,  for  small  sampling  times, 
the  simplified  performance  index  proved  to  provide  a  very 
accurate  approximation. 
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DESIGN  OF  A  DIGITAL  CONTROLLER 


FOR  A  LARGE  FLEXIBLE  SPACE 
STRUCTURE 

I  I ntroduc  ti  on 

The  promise  of  a  viable  space  transportation  system 
in  the  near  future  has  inspired  a  new  breed  of  satellite 
designs.  These  new  designs  are  far  larger  and  structurally 
more  flexible  than  their  predecessors.  Flexible  structures 
of  this  type  can  present  problems  to  the  control  designer  as 
they  may  exhibit  hundreds  of  low-frequency  vibrational  modes 
of  oscillation .  Using  classical  control  methods  it  is 
extremely  difficult  to  design  a  controller  to  insure  the 
stability  of  so  many  modes. 

One  of  the  most  widely  accepted  design  techniques 
to  counter  this  problem  is  the  application  of  modern  state 
space  control  theory  to  a  reduced  order  controller.  This 
technique  can  be  applied  to  a  wide  variety  of  flexible 
structures.  It  has  previously  been  applied  using  optimal 
time-invariant  linear  regulator  theory  for  active  control 
(Ref  1). 

The  active  control  is  applied  to  only  a  small  number 
of  critical  modes,  called  the  "controlled  modes".  The  choice 
of  which  modes  to  control  is  based  on  desired  system 


response.  Usually  the  lowest  frequency  modes  are  selected; 
however,  this  is  not  always  possible  nor  always  desired. 

While  the  controller  acts  on  the  controlled  modes,  the  con¬ 
trol  inputs  and  sensor  outputs  are  affected  by  the  remaining 
modes.  These  effects  are  known  as  "control  spillover"  and 
"observation  spillover",  respectively  (Ref  2).  The  effects 
of  some  of  these  uncontrolled  modes  can  make  the  system 
unstable.  Therefore,  the  effects  of  such  modes  must  be 
"suppressed"  in  the  controller  design.  These  modes  are 
known  as  "suppressed  modes".  The  remaining  modes  are  not 
modeled  in  the  controller  design  and  are  called  "residual 
modes" . 

In  previous  work  (.Ref- 1),  this  design  technique  has 
been  applied  in  the  continuous-time  domain.  However,  any 
control  system  based  on  a  digital  computer  must  be  a  discrete¬ 
time  system.  Meirovitch  (Ref  3)  has  shown  that  discretizing 
a  system  designed  in  continuous-time  can  lead  to  instabili¬ 
ties  as  the  sample  time  gets  large.  To  avoid  this  problem, 
such  a  system  should  be  designed  in  discrete  time. 

The  principle  purpose  of  this  thesis  is  to  design  a 
digital  controller  for  a  representative  space  structure. 

The  effects  of  the  suppressed  modes  and  of  different  sampling 
times  will  be  examined.  Two  variations  will  be  made  to  the 
design  of  the  controller  and  a  range  of  sampling  times  will 
be  used.  Since  many  satellites  are  used  to  ’.  earn  energy  or 
information  back  to  the  Earth,  the  performance  of  the 
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different  designs  will  be  judged  using,  as  a  figure  of 
merit,  the  deviation  of  a  nominal  line  of  sight. 


The  first  design  variation  involves  the  effects  of 
the  suppressed  modes.  The  system  will  first  be  designed 
ignoring  the  suppressed  modes.  The  system  will  then  be 
designed  to  eliminate  the  spillover  terms  in  the  state  space 
matrix  so  as  to  suppress  the  modes.  The  performance  of  the 
two  designs  will  then  be  compared. 

The  second  design  variation  involves  the  performance 
index  used  to  obtain  the  feedback  gains.  At  first  a  simpli¬ 
fied  approximation  to  the  discretized  performance  index  is 
used.  Then  a  much  more  accurate  performance  index  is  used 
and  the  two  systems  compared. ' 

The  first  sampling  time  examined  is  one  whose  sam¬ 
pling  frequency  is  well  above  that  of  the  highest  natural 
frequency  of  the  modeled  modes.  This  sampling  time  should 
provide  good  system  response.  A  range  of  larger  sampling 
times  will  then  be  used  to  see  the  effects  of  sampling  time 
on  system  stability. 


1 1  S  y s  L om  Model 


General  Description 

Since  a  mathematical  analysis  of  a  continuously  dis¬ 
tributed  space  structure  is  computationally  impossible,  the 
structure  is  modeled  as  a  discretely  distributed  system.  The 
model  used  is  one  developed  by  the  Charles  Stark  Draper 
Laboratory  (Ref  4).  It  was  developed  to  act  as  a  standard 
space  structure  for  controller  designs.  The  model  is  a 
tetrahedron  composed  of  lumped  masses  and  interconnected 
truss  members.  The  truss  members  are  assumed  to  be  massless 
and  exert  no  bending  moments,  only  axial  forces.  This  model 
is  illustrated  in  Figures  1  and  2.  It  consists  of  ten  nodes 
at  the  numbered  grid  points.  The  node  coordinates  are  given 
in  Table  I.  The  masses  are  located  at  nodes  one  through 
four.  Each  mass  has  three  degrees  of  freedom,  which  gives 
the  system  twelve  modes.  The  nodes  at  grid  points  five 
through  ten  are  fixed.  Six  collocated  sensor/actuator  pairs 
are  located  on  the  truss  members  extending  from  these  nodes 
and  are  oriented  along  the  members.  The  sensors  sense  posi¬ 
tion  only. 

The  key  results  of  a  NASTRAN  eigenvalue  analysis  of 
this  model  are  presented  in  Table  II.  The  eigenvectors 
associated  with  these  eigenvalues  are  shown  in  Appendix  A. 

For  the  time  response,  the  nominal  line  of  sight  of 
the  tetrahedron  shaped  "antenna"  is  along  the  z-axis  from 
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node  1.  Since  displacements  along  the  z-axis  do  not  dis¬ 
turb  this  line  of  sight,  only  perturbations  in  the  x-  and 
y-directions  are  considered.  The  mode  initial  conditions 
are  listed  in  Table  III. 


Table  I 

Node  Coordinates 


Node 

X 

Y 

Z 

1 

0.0 

0.0 

10.165 

2 

-5.0 

-2.887 

2.0 

3 

5.0 

-2.887 

2.0 

4 

0.0 

5.7735 

2.0 

5 

-6.0 

-1.1547 

0.0 

6 

-4.0 

-4.6188 

0.0 

7 

4.0 

-4.6188 

0.0 

8 

6.0 

-1.1547 

0.0 

9 

2.0 

5.7735 

0.0 

10 

-2.0 

5.7735 

0.0 

Mode 

Key  Results  of 

General ized 
Mass 

Table  II 

NASTRAN  Eigenvalue 

Generalized 

Stiffness 

Analysis 

wn( rad/sec) 

1 

1.0E+00 

1. 37E+00 

1. 17E+00 

2 

1.0E+00 

2.15E+00 

1.47E+00 

3 

1.0E+00 

8. 79E+00 

2.96E+00 

4 

1.0E+00 

1.26E+01 

3. 56E+00 

5 

1. 0E+00 

1.48E+01 

3. 85E+00 

6 

1.0E+00 

2.65E+01 

5.15E+00 

7 

1.0E+00 

3. 22E+01 

5.67E+00 

8 

1.0E+00 

3. 26E+01 

5. 71E+00 

9 

1. 0E+00 

7.99E+01 

8.93E+00 

10 

1.0E+00 

1.06E+02 

1 .03E+00 

11 

1.0E+00 

1. 19E+02 

1.09E+01 

12 

1.0E+00 

1 .95E+02 

1.40E+01 
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Table  III 


Initial  Conditions  for  Time  History  Response 


Mode 


)laccment 


Velocity 


1 

-0.001 

-0.003 

2 

0.006 

0.010 

3 

0.001 

0.030 

4 

-0.009 

-0.020 

5 

0.008 

0.020 

6 

-0.001 

-0.020 

7 

-0.002 

-0.003 

8 

0.002 

0.004 

9 

0.000 

0.000 

10 

0.000 

0.GC0 

11 

0.000 

0.000 

12 

0.000 

0.000 

Equations  of  Motion 

The  vibrational  motion  of  a  large  space  structure 
is  given  by  the  equation 


C\  |U)  +■  E  ^(t)  +  K  lit)  =  Bii«) 

where 

q  =  an  n-vector  of  generalized  coordinates 
M  =  an  nxn  symmetric  mass  matrix 
E  =  an  nxn  symmetric  damping  matrix 
K  =  an  nxn  symmetric  stiffness  matrix 


u 

D 


an  nac^.-vector  of  inputs 

an  nxnac^_  matrix  of  actuator  coefficients 


and  where  the  dimensions  n  and  n„  ,  are  the  number  of  modes 

(XC  c 


and  the  number  of  actuators,  respecitvely . 


Equation  (1)  can  be  written  in  modal  coordinates,  n 


by  using  the  transformation 
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<£(*)  -  + 

where  <p  is  the  nxn  modal  matrix  for  Eq  (1).  This  modal 
matrix  satisfies  the  conditions 

<)>TAl  j>  -  I 
4>T  K  <j>  = 
cj>T  E  <J>  —  2  "5 

wli'ere 


CO 


(3) 

(O 

C.O 


I  =  an  nxn  identity  matrix 
2£<j j  =  an  nxn  diagonal  damping  matrix 

p 

co  =  an  nxn  diagonal  matrix  of  eigenvalues  of  Eq  (1) 
With  this  transformation,  Eq  Cl)  becomes 


not)  +-  vo  +  w*  ^t)  -  co 


Equation  (6)  put  into  state  vector  form  is 

X  Ct)  -  A  x(t)  t  ft  u(t)  (7) 


where 


A  = 


-oj1  ["-2  3^ 


co 

n) 

(to) 
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The  output  equation  from  the  sensors  is  of  the  form 


yCO  =  Cf  \(.t)  <") 

where  C  is  an  n  xn  matrix  of  sensor  coefficients .  The 
p  sen 

dimension  n  is  equal  to  the  number  of  sensors, 
sen  1 

Using  the  transformation  of  Eq  (2),  the  output 
equation  becomes 

y  tt)  =  C  xU)  <l2> 

where 

C  =  [cr+  ;  o] 

Equations  (7)  and  (.12)  are  the  equations  of  motion 
in  the  form  to  be  used  in  the  controller  design.  The  state 
vector,  x,  is  of  dimension  2n .  As  previously  discussed, 
the  system  modes  can  be  separated  into  controlled,  suppressed, 
and  residual  modes.  When  such  a  partitioning  is  made,  the 
corresponding  states  of  the  state  vector  are  also  partitioned. 
The  controlled  states  are  those  states  directly  controlled 
to  achieve  the  specified  performance.  The  suppressed  states 
are  those  uncontrolled  states  which  may  still  strongly 
influence  the  system’s  performance  and  should  be  accounted 
for  in  the  design  of  the  controller.  The  residual  states 
correspond  to  those  modes  which  have  such  a  high  natural 
frequency  that  they  will  probably  not  be  excited  by  the 
controller  and  thus  may  be  left  out  of  the  controller  design. 
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With  this  partitioning  of  states,  the  state  equation  becomes 


where 


icct)  =  +  6C  m  0*M 

xsct)  =  As  xsa)  +  6,  /tC*)  (ml) 

xr(i)  -  Ar  xrU)  +-  lbrMtt) 


i 


r 

i 


x 


C»SA 


:  C,S 


i  I 


(/st) 


are 


B 


c 


T 

The  matrices  <(>  .  D  are  those  matrices  whose  elements 
J 

linear  combinations  of  the  mode  shapes,  <f> .  The  matrices 
Bs>  Br  are  shown  in  Appendix  B, 

With  partitioning  the  output  equation,  Eq  (12) 


becomes 


y(*'>  =  Cc  5^*1  +  c.sx,ut)  +crzr(t) 


(ii) 


where 


c<  -  cp  4>-  1 

j  L  p  ’J  • 


o 


J - 


Cio) 


The  matrices  C  <f>  .  are  also  matrices  whose  elements 
P  J 

are  linear  combinations  of  the  mode  shapes.  For  the  special  case 


11 


whore  tho  position  sensors  and  point  actuators  arc  collo¬ 
cated  and  aligned,  the  B.  and  C.  matrices  are  related  by 

J  J 

c  =  rt  t”-) 


Cs  =  tJn> 

C  -  &'r  (21c.) 

•  r 


In  order  to  design  a  digital  controller  in  discrete 
time,  these  equations  must  be  descritized. 

Equations  of  Motion  -  Discrete  Time 

In  a  sampled-data  control  system  the  control  inputs, 
u^,  instead  of  changing  continuously  in  time,  are  incre¬ 
mented  only  at  discrete  time  instances.  Between  those 
instances,  the  input  is  held  constant.  Thus,  the  discrete 
control  input,  u^,  is 

U  (*)  =  -u.(kT)  kT  <  jt-  <  (kii)T  C?-0 

where  k  =  0,  1,  2,  ...,  and  T  is  the  sampling  time. 

It  can  be  shown  that  the  state  space  representation 
of  the  system  in  discrete  time  is  then 

x  E(V*  i)Tl  -  £(T)  x(\<T)  +  PCt)  /a (\rr)  ( ?h) 


12 


Tho  matrix  i>(T)  is  related  to  the  continuous-time  A  matrix 
by 


I(T)=  e"T 


(2S) 


F(T)  is  related  to  the  continuous- time  B  matrix  by 

r(T)  =  |re,T  jr  B  C3 4) 

/<i 


$(T)  and  T (T )  are  not  functions  of  time  but  are 
functions  of  the  system  sampling  time.  This  sampling  time 
is  usually  a  constant  in  the  system  as  it  will  be  in  this 
analysis.  Therefore,  HT)  and  T(T)  are  constants  and  for 
convenience  will  be  written  as  $  and  T. 

As  in  the  continuous-time  domain,  the  state  equation, 
Eq  (24),  can  be  partitioned  into 


ScftW  OTJ 

x.Hu+rt-rJ 

X,[UvOt] 


where 


r 

c 


Ic*cdkT)  +  Pt  u  (.VT) 

xs(kT)  4-  rs  m  CkT) 


Ir  xr  CWT)  -b  £  JU  CkT) 


c. 

A.r 


r '  A~ 

f *  e  r 


(.27  0 

(27  c) 


0*0 
t?  %  O 
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r 


1..=- 

eA'T 

Ol«) 

n  = 

CT  A.T' 

XT 

(2  11) 

4  - 

eA-T 

(l  oa) 

r 

E  = 

'  o 

Jr 

K> 

r 

(3  ol) 

The  continuous-time  output  equation,  Eq  (19), 

becomes 

y  (  kT)  =  CUT)  +  cs  As  CUT)  +-  Cr  Kr  0<T)  (31) 


Equations  (27)  and  (31)  are  the  equations  to  which 
the  controller  design  will  be  applied. 


In  all  of  the  following  designs,  the  feedback  gains 
for  the  controller  and  the  observer  are  designed  using 
optimal  control  theory.  There  are  some  differences  between 
the  performance  index  for  a  discrete-time  system  and  that  for 
a  continuous-time  system.  In  this  chapter,  discrete  optimal 
control  theory  is  reviewed  to  demonstrate  these  differences. 

In  the  cont inuous- time  case,  the  performance  index 
frequently  used  for  systems  described  by 

X  (X)  —  A  X  C t )  -h  G>  AX  Ct ) 

is  given  by 

J  i  5  55  x(t)  t  R,  utt)  At  (33j 

O  — 

where  Q  is  a  2nx2n  positive  semi-def inite  state  weighting 
matrix  and  R  is  an  n  .xn  .  positive  definite  control 

^CX  clCX 

weighting  matrix.  The  solution  to  this  problem  is  known 
to  be 

(X)  =  Fx(>)  (.3*1) 

where  F,  the  feedback  gain  matrix,  is 


F  —  -  K  V  P 


(35) 


and  where  P  is  the  solution  to  the  steady  state  Ricatti 
equation 


?^  i-  /Tp-  +-  ^  =  ° 


Czi) 


For  the  discrete-time  case,  an  appropriate  perfor¬ 
mance  index  for  systems  described  by 

*Lu*iyri  -  i.  xOro  i-  n>u,ci<T)  (37) 

is  given  by 

^  ~  *  2.  j  xTa)  Gl  x  Li)  +  (})  KmCX)  dt  (3O 

'  kT  - 

But  for  discrete  systems 

AbU)  =  <UC\<T) 

AU-WT)  /t  .  /  AC*-r) 

X(*)  =  €  X  CkT)  -v  (  e  (,ir  B  tj CkT)  ^ ) 

k“r 

for  kT  <  t  <  (k+l)T.  As  shown  in  the  preceding  chapter, 

Eq  (39)  can  be  written 

xa)  -  \Lf-vj)  xtkr)  +  ra-kT)  xxOcT)  (40^ 

After  making  the  transformation  r=t-kT,  the  performance 
index  becomes 


16 


I  00  f 

J  - ^  Jo  [£'rl  2  (i(T)  +  [}(T)  x CP)  + 

rtr)  +  2To<-t)  jr 

)  00  7" 

=  1  (  C*^t)  lx(r)J  q  JVrj  X  ()<T)]  +- 

[xr0<T;  q  £  r(T)  m  ckT)]  + 

tT  d^kT)  PT(r)J  Q  [  f  (r)  <a<T)J-h 

L  -^T^kT)  PT(t)3  <3  Cn  (T)  A4  (.kT)]+- 

M.TCkT)  K  xaOcT)  AT  Cm) 

Finally,  the  performance  index  can  be  written 

T  =  i  f  xT0<r)  Sj  x,  CkT)  +  xTCkT)  sj  MCkT)  + 

XlTCkT)  Sj~  X  (kT)  4-  MT CkT)  RjX/CkT)  C<I2) 

where 

^  =  $  1T(*)  Q  &(*) 

o 

sa  =  ^  iTa)  Q  net)  JJr  O) 

=  $T  rTct )  <s*  na)  +-  r  (Mr) 

Equations  (43)  through  (45)  can  be  simplified  by  the 
approximation 
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AX  _  t  j_  a  +■  l  A7  X  A"3 X 

c  -  1  +  -57-  ■»-  -TT  +  ,  .  . 


It  can  be  shown  that,  if  the  product  of  T  with  the  largest 
eigenvalue  of  A  is  very  small,  that  is  if 

|  _/V.  T  |  <<  1 

1  mox  » 


where  A  is  the  largest  eigenvalue  of  A,  then  Eq  (46) 

ulciX 

can  be  approximated  by 

AX  _ 

e  -  -L 

With  this  approximation,  Eqs  (43)  through  (45)  become 

-  sr 

Sj  =  O 

and  the  simplified  performance  index  can  be  written 

J  -  i  I  X7- (kT)  QT  X  (Vt)  -f  xT| KT  M  (kT) 


For  a  state  equation  of  the  form  of  Eq  (37),  the 
solution  for  this  simplified  performance  index  is 


(*n) 


(4  g) 


(HI) 

(50) 

(51) 


(.52) 


^CkT)  -  F  x  CkT) 


whore 


F  =  -  ]T  1  +  Crt)~'  rT  f  r  J  '(  rt )-1  p1-  p 


(E3) 


(54) 
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and  where  P  is  the  solution  to  the  discrete  steady-state 
Ricatti  equation 

^  ~  -  1T  pp  (rt  t-  r^rry1  pxp  $  +  c?t  css) 
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IV  Design  _I 


System  Design 

The  first  design  makes  the  assumption  that  the  sup¬ 
pressed  modes  can  be  ignored.  The  simplified  approximation 
of  the  performance  index  is  used.  The  state  equation  for 


the  controller  is 

xJLcv+n-r]  =  \kcc\cT)  +-  n  m  ci<t)  (56) 

and  the  output  equation  is  assumed  to  be 

y  CkV)  =  xoCWT)  (57) 

The  performance  index  for  this  system  is 

Oc=> 

*I(VT)  <$.T  XctW  T)  f  MT  (kT)  R.T  WkJ)  CsO 

k=i 

From  Chapter  III,  the  solution  to  this  problem  is 

AX*(XT)  -  F  XtCkT)  C51) 

where 

F  —  ~  Ur  +•  CrtF'  rrfpj  'crt)-1  pcrp  ^  (to) 


and  where  P  is  the  solution  to  the  discrete  steady-state 
Ricatti  equation 


20 


P  =  l J"  p  Ckt  +-  nTrr  )  '  f^Tp  5^  +  ax  0 

Then  the  state  equation  can  be  rewritten 

^rCk+0T]=5  X,tkT)  *-  r  F  *  CkT)  (<5.2 ) 

*  C*  C»  *■  C  C. 


The  states  are  not  directly  available  to  the  con¬ 
troller;  therefore,  an  observer  is  required  to  provide  in¬ 
puts  to  the  controller.  The  observer  state  equation  is 

j^Da.OT]  =  I^OcT)  t  rc  <d(Xr)  i-  KCjOtV-  yCkT)]  Ui) 

where  K  is  the  observer  feedback  gain  matrix  and  the  observer 
output  equation  is 

y  Cut)  -  Cc  xcO<t)  cch) 

The  output,  y(kT),  is  still  assumed  to  be  determined  by 
Eq  (57).  Thus,  Eq  (63)  can  be  written 

>LE(u.ot1=  C^-kcci  qF)xt(kT)  +  KCckcO<V  US) 

To  determine  K,  a  new  state  vector  must  first  be 
defined;  the  error  state  vector 

eUT)  H  X  CkT)  -  X  O ;T)  LU) 

Then 
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ejCCk+DT 7  =  XttUtlYT]  -  xJUvOTj 

=  x^lj) KC^i  P  f]  XtO<X) 

-  |<C  xca<-r) 

—  c  50-KCc)  e  0<T)  (67) 

It  is  desired  to  design  K  such  that  e  is  minimized. 

Eq  (67)  is  not  of  the  form  of  Eq  (37)  from  Chapter  III,  for 
which  the  optimal  solution  is  known.  However,  the  eigen¬ 
values  of  a  matrix  are  the  same  as  those  of  its  transpose. 
Therefore,  the  following  equation  will  have  the  same  closed- 
loop  eigenvalues  as  Eq  (67). 

w[Cic*OT]=  £>0<t)  -  c^}CkT)  Ui*) 

where 

(_\<7)  -  KT^OcT)  U*\) 


This  equation  is  of  the  form  of  Eq  (37)  so  that  the  perfor¬ 
mance  index  is 

wT(kT)  ^  T  y^T)  t  j.TCU T)  RoT  (7c>) 

The  solution  is 

«*0<r)  =  X^uCkT) 

o 
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whore 


KT=  [it-  CR.T)''ct 


?  c 

°e>  c 


■J' 


p  I 


°E 


(7-?) 


and  where  PqQ  is 
Ricatti  equation 


the  solution 


to  the  discrete  steady-state 


P  -I 

OR,  « 


1  -  ? 


r  c 


OL/r-i- 


C  f  c 

c  1  op,  * 


*'  +  <3' 


(73) 


T 

The  transpose  of  K  gives  the  desired  observer  feedback  gain 
matrix,  K. 

Thus  the  system  equations  based  on  the  assumption 
that  the  suppressed  modes  can  be  ignored  are 

X^DuOTl  -  XcCl<T)+  PF  xJV T) 


—  P  f)  )  —  )~^  j- e  Ck"1")  (74) 

e  rCk*0T ]  =1  Uc'KCt)e(^)  (75) 

In  many  cases,  the  suppressed  modes  cannot  be  simply 
ignored.  To  see  their  effect,  the  total  closed  loop  system 
of  equations  is  derived.  The  output  equation  was  assumed 
to  be  Eq  (57);  however,  with  the  suppressed  modes  included, 
it  is  actually  (ignoring  the  residual  modes) 

y(kT)-  CcxtCkT)+  Cs  xsOcT) 

With  this  output  equation,  the  observer  state  equation 
becomes 
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Aja.OT]  =  (i<-  KCtt  riF)£ar)i  KC^xJk-r)  +  KO^CVT)  07) 

The  error  state  equation  is  thus 

c[ Jk+Orl  -  C£c~l<0  e(kT)-  KC,  ^0<T)  C-n) 

and  the  suppressed  state  equation  is 
*sfrktl)T]  =  £5  x/)<T^  -t-  f] 

=  xsCk~r)  t-  rs  xcO<T) 

-  QPx/kT)-  Ps  F  e0<0  +  £sX/kT0  (71) 

The  total  system  state  equation  can  be  written  in 


matrix  form  as 


xJLOc+DT] 

j  tpcF  1  -  r  f  1  0 

c  c  1  c  1 

■  XJU)" 

c  [CkiorJ 

— 

1  1 

O  1  1  -  kcc  1  kc5 

1  t 

e  CkT) 

x^rtwnrj 

n  f  1  -  r.  f  1  ? 

L  .  1  ‘  1  s  J 

*say) 

The  T  F  and  KC  terms  are  the  control  and  observa- 
s  s 

tion  spillover  terms.  The  eigenvalues  of  !I>C  +  TcF  , 

$  -  KC  ,  and  $  are  the  designed  eigenvalues.  Unfor- 
tunately,  due  to  the  spillover  terms,  the  total  system 
eigenvalues  are  altered  from  the  desired  eigenvalues.  As 
the  gain  matrices  F  and  K  are  increased  (with  increased 
control  of  the  system),  the  altered  eigenvalues  can  become 
unstable . 
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Sys  tein  Per  romance 

In  the  following  analysis,  the  system  eigenvalues  are 
plotted  in  the  z-domain.  Before  examining  the  system  eigen¬ 
values  and  time  response,  it  is  perhaps  best  to  review  the 
interpretation  of  eigenvalue  placement  in  the  z-plane.  There 
is  a  correlation  between  eigenvalue  position  in  the  s-plane 

and  that  in  the  z-plane.  The  s-plane  is  mapped  into  the 

sT  aT 

z-plane  via  the  transformation  z  =  e  =  e  L  ±  w^T  where 
s  =  a  ±  .  This  transformation  maps  the  entire  left- 

half  s-plane  into  a  unit  circle  centered  at  the  origin  in 
the  z-plane.  In  order  for  a  system  to  be  stable,  all  eigen¬ 
values  must  lie  within  the  unit  circle.  The  entire  real 
axis  in  the  s-plane  is  mapped- onto  the  positive  real  axis 
in  the  z-plane.  Lines  of  constant  damping,  a,  in  the  s-plane 
map  into  concentric  circles  centered  at  the  origin.  a=0 
lies  on  the  unit  circle  and  a=-»  is  at  the  origin.  Lines 
of  constant  frequency,  map  into  radial  lines.  And  lines 

of  constant  damping  ratio  become  logarithmic  spirals  in  the 
z-plane . 

Listed  in  Table  IV  are  the  eigenvalues  of  the  above 
system  design  for  a  sampling  time  of  0.1  seconds  and  with 
the  state  and  control  weighting  matrices  entered  as  identity 
matrices.  The  system  eigenvalues  are  plotted  in  the  z-plane 
in  Figure  3.  The  time  response  to  the  initial  conditions 
given  in  Chapter  III  is  shown  in  Figures  4  and  5. 
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As  a  basis  Cor  comparison,  the  discrete  open- loop 
eigenvalues  and  time  response  are  presented  in  Appendix  C. 

It  can  be  seen  that  the  open- loop  system  is  stable  but  the 
transient  response  is  unsatisfactory.  After  twenty  seconds, 
the  X  and  Y  displacements  of  the  line  of  sight  are  still 
quite  large.  For  this  design  and  those  which  follow,  it  is 
desired  to  reduce  the  X  and  Y  displacements  to  less  than 
0.0004  radians  and  0.00025  radians  respectively,  in  twenty 
seconds . 

Table  IV  demonstrates  the  effects  of  the  spillover 
terms.  The  designed  eigenvalues  are  listed  in  the  left-hand 
column  and  the  resulting  system  eigenvalues  from  Eq  (80) 
are  listed  on  the  right.  Because  the  weighting  matrices  Q 
and  QqB  are  only  identity  matrices,  the  gain  matrices  F  and 
K  are  small  and  the  system  eigenvalues  are  not  altered 
drastically.  The  system  is  still  stable  and  Figures  4  and 
5  show  that  the  time  response  has  been  improved.  The 
desired  specifications  are  still  not  met,  however. 

To  meet  the  specifications,  more  emphasis  must  be 
placed  on  controlling  the  states.  This  is  done  by  increas¬ 
ing  the  state  weighting  matrices  Q  and  in  the  perfor¬ 

mance  indices.  Table  V  presents  the  eigenvalues  of  the 
system  design  for  Q  and  increased  to  10  I  matrices. 

For  this  case,  the  spillover  terms  are  increased  enough  to 
drastically  change  the  designed  eigenvalues.  In  fact,  one 
of  the  suppressed  eigenvalues  is  now  unstable.  These 


system  eigenvalues  are  also  plotted  in  Figure  3.  The  time 
response,  shown  in  Figures  6  and  7,  bears  out  this  insta¬ 
bility.  Further  trials  of  Q  and  would  show  that  this 

controller  design  cannot  achieve  the  desired  system  perfor¬ 
mance  because  the  system  goes  unstable  before  the  desired 
performance  can  be  reached. 
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Table  IV 


Eigenvalues 

of  Design 

I 

Q  = 

1.0  I  ;  q0B  . 

=  1.0  I;  T 

=  .1 

Eigenvalues  of 

+  roF 

System  I 

^  i  <r 

onvalues 

.98548  ±  . 

11585i 

-C- 

.98572 

f 

. 11690i 

.91600  ±  . 

33968i 

-c- 

.91587 

+ 

. 33988i 

.97811  ±  . 

14434i 

-c- 

.97810 

+ 

. 14 54 7 i 

.90144  ±  . 

36415i 

-c- 

.90142 

+ 

.  364 16 i 

Eigenvalues  of 

$  -  KC 

c  c 

.92264  ±  . 

10453i 

-0- 

.92237 

+ 

. 09049i 

.76201  ±  . 

22371i 

-0- 

.76207 

+ 

.223001 

.88784  ±  . 

111419i 

-0- 

. 88937 

+ 

. lOllli 

.72520  ±  . 

21495i 

-0- 

.725 40 

± 

.214701 

Eigenvalues  of 

.95496  ±  . 

29170i 

-s- 

.95363 

± 

.  29292i 

.84081  ±  . 

53607i 

-s- 

.84090 

± 

. 53862i 

.86809  ±  . 

49121i 

-s- 

. 86809 

± 

.491211 

.83892  i  . 

53899i 

-s- 

. 83856 

± 

. 54050i 
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LOSY  *10‘" 
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Table  V 


Eigenvalues  of  Design  I 
Q  =  10.0  I  ;  Qqb  =  10.0  I  ;  T  =  .1 


Eigenvalues 

of  $  +  r  F 

c  c 

System  3 

Eigenvalues 

.96930  ± 

. 11360i 

-C- 

.97633 

+ 

. 12079i 

.87318  ± 

. 31781i 

-C- 

.87073 

+ 

.  31938i 

.95475  ± 

. 13946i 

-C- 

.95827 

+ 

. 15080i 

.85127  + 

.  33561i 

-C- 

.85064 

+ 

.335751 

Eigenvalues  < 

of  $  -  KC„ 

c  c 

.29844  + 

Oi 

-0- 

.29965 

+ 

Oi 

.87953  + 

Oi 

-0- 

.85149 

+ 

Oi 

.80409  ± 

.  01949i 

« 

-0- 

.98221 

+ 

Oi 

.23649  + 

Oi 

-0- 

.23670 

+ 

Oi 

.58986  + 

Oi 

-0- 

.61812 

+ 

Oi 

.88935  + 

Oi 

-0- 

.84875 

+ 

Oi 

.84675  + 

Oi 

-0- 

.53363 

+ 

Oi 

-0- 

.95642 

+ 

Oi 

Eigenvalues  < 

of  $s 

.95496  ± 

. 29170i 

-s- 

.94752 

± 

.29673i 

.84081  ± 

.  53607 i 

-s- 

.85356 

+ 

.  57269i 

.86809  ± 

.  49121i 

-s- 

.86809 

+ 

.  49121i 

.83892  ± 

.53899 

-s- 

.83686 

+ 

.  55492i 
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Figure  6.  I, CSX  VS.  TIFE,  Design  I,  Q  =  Qr,p  =  10  I,  T  =  .1 
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0.06 


V  Design  II 


In  the  second  design,  the  effects  of  the  suppressed 
modes  are  accounted  for.  The  system  is  designed  to  elimi¬ 
nate  their  influence,  hence  their  effects  are  "suppressed". 


The  design  begins  with  the  same  state  equation  as  the  first 
design 

but  this  time  the  output  equation  is  assumed  to  be 

y  CIO-)  -  C,_  Xc(.lcT)  +-  cs  x^CkT)  (sO 

The  performance  index  is  the  same  as  before 

_ 

I  =  i  Z.  *>t)  x/kT)  +-  /Ux<:kT)  kTxxC\cT)  cn) 

k-( 

The  solution  is  again 

,U*~(kT)  -  F  X^CkT)  (H) 

This  time,  however,  in  order  to  eliminate  the  control  spill¬ 
over,  F  must  satisfy  the  following  conditions 

rt  f  *  o  f  *=) 

Hs  F  -  O  Lit,) 


This  can  be  achieved  by  defining  a  new  control  variable,  v 
linearly  related  to  u  by 

-  n  vcwt) 


35 


where  the  matrix  M  satisfies  the  relationships 

nc  n  *  o  (s^ 

r3  h  -  o  c*i) 

By  applying  optimal  control  theory  using  v  as  the  control 
vector,  a  control  law  of  the  form 

V(k~r^  =  6  XcO<T)  o ) 

can  be  found.  Transforming  back  to  u  one  finds 

XL  CkT)  —  /16*eUT)  =  FxeO( T) 

which  satisfies  Eqs  (85)  and  (86).  Applying  the  transforma¬ 
tion  given  in  Eq  (87),  the  state  equation  becomes 

Xer^H)T3  =  lcxtCkTj  t  rcM  v  Cut )  no 

and  the  performance  index  becomes 

T=  i  f.  xJCkT)  QT  X  t(kT.  VT0<T)  ^  v  UcT)  (10 

K -I 

T 

where  R^  =  M  RTM  .  The  solution  to  this  problem  is 

Y*OT)  =  G  x.XW'O  (13) 

where 

G  =  -  ti+-  R>TncTpn  a.]*1  rT  p  Ci‘0 

and  P  is  the  solution  to  the  discrete  Ricatti  equation 
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p  -  ~  M  Ckm  +■  //Tr'rpr  /i  )~]/iTrc'r r  *- 

The  matrix  M,  which  satisfies  the  conditions  of  Eqs 

(88)  and  (89),  is  found  by  the  method  known  as  Singular 

Value  Decomposition  (Ref  4).  It  is  desired  to  satisfy 

Eq  (88)  with  a  non-zero  matrix  M.  Since  T  is  given  by 

s 

Eq  (29b),  Eq  (88)  can  be  written 

,  kr  A  r 

O 

Thus,  to  satisfy  Eq  (86)  it  is  only  necessary  that 

B  srt  =  o 

In  Chapter  II,  it  was  shown  that  Bg  is 


where 


Thus,  to  satisfy  Eq  (97),  it  is  only  necessary  that 


e.  n 


o 


Now,  using  a  singular  value  decomposition  of  the  matrix  B 

~  S 

the  matrix  can  be  written  as 


B>  - 


3T  (IS) 


0~0 


no 


(01 ) 

(loo) 


(>0  \) 


U  HY 


where 


U  is  an  nsxng  orthogonal  matrix  of  left  singular  vec¬ 
tors 

V  is  an  n  ,xn  orthogonal  matrix  of  right  singu- 
act  act 

lar  vectors 

and 

N  | 

o  \  o 


such  that  N  is  an  n^xn^  diagonal  matrix  of  the  singular 


values  of  B  .  The  dimension  n„  is  the  number  of  non-zero 
~s  a 


singular  values.  Each  singular  value  is  real  and  must  be 
greater  than  or  equal  to  zero:  The  total  number  of  non-zero 


singular  values  is  equal  to  the  rank  of  B  .  If  B  is  full 

~  o  ~  o 


rank,  then  nQ  is  equal  to  either  ng  or  nact  ,  whichever  is 
least.  When  the  number  of  actuators  is  greater  than  the 
number  of  suppressed  modes  (which  is  the  case  for  this 


analysis),  then  n  <  n  ,  and  n  =  n  .  In  this  case,  the 

s  act  o  s 


matrix  E  becomes 


2  =  [n  j  °  ] 


Now,  the  matrix  V  can  be  partitioned  into 


v  -  [y  !  vj 


where  Vi  is  n  ,xn  and  V„  is  n  ,x(n  ,-n  ).  Then  Eq  (100) 
x  act  s  «  act  act  b 


can  be  written 


or 


u 

'n  ;  o] 

V1' 

1 

1 

VT 

L.  *  J 

(loS) 


=  UNV|T  (io&) 

If  both  sides  of  Eq  (106)  are  multiplied  by  V2  ,  then 

=  trNY,TV, 

-  o  (10 -I) 

because  V  is  an  orthogonal  matrix.  The  matrix  V2  meets  the 
condition  of  Eq  (99),  thus 

M  =  V*  C)ot) 

is  a  possible  choice. 

The  state  equation  can  now  be  written 

xjXwo-O  -  +  rtp  xjat)  ooO 

where  F  is  given  by  Eqs  (90b),  (94),  and  (108).  The 
observer  state  vector  in  the  above  equation  is  given  by  the 
observer  state  equation 

xJUiOTl  =  +  rc  v  ytwr)  "  yCk"^  CMo) 

where  K  is  the  observer  feedback  gain  matrix  and  the  observer 
output  equation  is 

/(kT)  -  Cc  Xt(]<T)  OH) 
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Substituting  Eqs  (82)  and  (111)  into  Eq  (110),  the  observer 
state  equation  can  be  written 

?cr0c+l)Tj  =Ci^KCt  +  rtF)  AC0<T)  t  )<CoXc(),T)+  )<CS  X/)<T)  (1/2) 

To  determine  K,  the  error  state  vector  is  defined 

eCWT)=  /cCUT)  -  Xc  OcT) 

Then 

e[Cktl)T]  -  -  X.P^iOT] 

—  le  d/-1'1)  +  rc  F  Xc(kl)  -  (  tP  F)  Xt0<T ) 

-  KCc  Jt0(T)  -  KC,  x/WT) 

=  (ic-Kcci  eClcT)  -  |CC5  x/ kl)  CM 3) 

In  order  to  eliminate  the  observation  spillover,  the  gain 
matrix,  K,  must  satisfy  the  conditions 

O 

Kcs-o 

As  with  F,  the  observer  gain  matrix  can  be  defined  as  the 
product  of  two  matrices 

K  =  L/1oa 


(IM) 

015) 

( //(-) 
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whore 


C 1 1  7  ^ 


*«>,  cc  *  o 

^  =  o  (''O 

Then  the  error  state  equation  becomes 

e  runiTl  =  (?c-kcJ  eOcT)  On) 

As  in  Design  I,  this  equation  will  have  the  same  closed-loop 
eigenvalues  as  the  equation 

wltk+OT]  =  ?TdO<T)  -  CjjCkT)  (l?°) 

where 

jCkT)^  0<t;  tm) 

With  K  defined  as  in  Eq  (116),  this  equation  can  be  trans¬ 
formed  into 

wlA+OT]  -  iJwO'r)-  KOtT)  (122) 

where 

Vl(kT)  -  LT  uCkT)  (  7  i) 


The  performance  index  for  this  state  equation  is 

J=  i  X.  tfCk-O  S  T  wtlcT)  +■  VlUT)  R  yCkT) 

k^\  0  -  -  v  - 

where 


K  T  A 

orj  ->C. 


( nr) 
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The  solution  is 


K*o<t')  -  ltu  (m 


U-;} 


where 


L1  -  T~t  +-  R  '  An.  C  ?  CT/1  T  T1  R  m  c 

C£>n  C  C  oa  -1  'yMJ  '  'op  , 


onu  H  c  ot> 

A1 


1 


C  127) 


and  where  P^g  is  the  solution  to  the  discrete  steady-state 
Ricatti  equation 


»6 


( l?s 


The  transpose  of  L1  given  in  Eq  (127)  is  the  desired  matrix 
L  in  Eq  (116).  M^g  must  be  determined  by  another  singular 
value  decomposition.  This  time  it  is  desired  to  satisfy 
Eqs  (117)  and  (118).  In  Chapter  II,  it  was  shown  that  Cg  is 


c, =  0*  i  °] 


where 


(  12  l) 


c  —  c 

>-  9  ^-p  T s 


(  13  o) 


To  satisfy  Eq  (118),  it  is  only  necessary  that 


/\.  C  =  0 

°£  — S 


(>31) 


A  singular  value  decomposition  of  the  matrix  C  gives 

~  S 


£=  u.52,X. 


where 


(  137) 


UQB  is  an  nsenxnsen  orthogonal  matrix  of  left 
singular  vectors 
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and 


V nn  is  an  n  xn  orthogonal  matrix  of  right  singular 
Uij  is  & 

vce  tors 


N  n 

1  ~ 

I  Cl 

— 

o  13 

i 

0 

1  0 

i  J 

(  /33) 


such  that  NqQ  is  an  nQxna  diagonal  matrix  of  non-zero  singu¬ 
lar  values  of  C  .  The  dimension  n  ,  the  number  of  non- 
—  s  a 

zero  singular  values,  is  equal  to  the  rank  of  C  .  If  Cg  is 
full  rank,  then  n^  is  equal  to  cither  ng  or  nf,Gn  ,  whichever 
is  least.  When  the  number  of  sensors  is  greater  than  the 
number  of  suppressed  modes  (which  is  the  case  for  this 
analysis),  then  n  <  n  and  n^  =  n,  .  In  this  case,  the 
matrix  becomes 


C 


03  «  ) 


Now,  the  matrix  UQB  can  be  partitioned  into 


u 


on 


U  1  u 

VOEI  |  O  liX 


C  bO 


where  UQB1  is  ngenxns  and  UQB2  is  nsenx(nsen-ns) .  Then 


Eq  (132)  can  be  written 


C  - 


6  \  1  o  ■>> 


>] 


M 


•v 


(  /36  ) 


Since  the  matrix  UqB2  is  multiplied  by  zero,  this  equation 
becomes 


4  3 


037) 


r 


£.=  u..,wv#s 


If  both  sides  of  Eq  (137)  are  multiplied  by  U 

XJT  C  -  UTV  N  "V 

<,&1.  ^  i  *51  Otfl  <05,  (3d, 


-  o 


0B2 


then 


C  13  f  ) 


T 

because  U^g  is  an  orthogonal  matrix.  The  matrix  U^gg  meets 
the  condition  of  Eq  (131),  thus 


yM  -  V 

'  oii  oe, * 


(m) 


The  error  state  equation  is  now  as  given  in  Eq  (119) 
with  the  observer  gain  matrix,  K,  as  given  by  Eqs  (116), 
(127),  and  (139). 


The  suppressed  state  equation  is 

aijcktorj  ~  J^CkT)  *  rs/jo<T) 

-  4S  XsO*0  rs  j=  Ac  ()<T ) 

=  i5Ks(kT) 


The  total  system  state  equation  can  be  written  in 
matrix  form  as 


xrEck+t  )t3 

%  n  p.  j  -  p  p  *o 

c  c.  1  1  c  I 

"xcoa) 

e  COctOT] 

— 

1  1 
o  i  ^cc  !  ° 

PUT) 

XjXk*  OtI 

r 

0 

0 

o 

i _ 

X  (U) 

(m) 
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System  Performance 


Table  VI  presents  the  closed-loop  system  eigenvalues 
for  Design  II.  With  the  spillover  terms  eliminated,  the 
system  eigenvalues  are  identical  to  the  desired  eigenvalues 
(to  at  least  eight  decimal  places).  The  eigenvalues  for 
Q  =  1.0  I  and  =  1.0  I  are  shown  in  the  left-hand 

column.  These  eigenvalues  are  plotted  in  Figure  8.  They 
are  all  stable  and  the  controlled  state  eigenvalues  (those 
of  4>  +  T  F  )  are  slightly  closer  to  the  origin  than  the 

v  C 

open-loop  eigenvalues  in  Appendix  C.  This  implies  that  the 
time  response  should  be  slightly  improved,  which,  as  seen 
in  Figures  9  and  10,  is  the  case.  In  comparison  with  the 
corresponding  system  response -of  Design  I,  it  can  be  seen 
that  the  eigenvalues  of  Design  II  are  slightly  farther  from 
the  origin  than  those  in  Design  I.  The  time  response  is 
actually  somewhat  poorer  for  Design  II.  Thus,  there  seems 
to  be  some  loss  of  control  (for  the  same  Q  and  QqB)  in 
return  for  suppression  of  the  spillover.  This  is  due  to  the 
transformation  of  Eq  (87),  because  the  new  control  vector, 


v,  is  of  lower  dimension  than  u  . 


To  improve  the  performance,  the  state  weighting 
matrices  Q  and  QqB  are  increased  by  a  factor  of  10.  The 
resulting  eigenvalues  are  shown  in  the  center  column  of 
Table  VI  and  are  plotted  in  Figure  8.  These  have  remained 
stable  (unlike  those  of  Design  I)  and  have  moved  closer  to 


the  origin.  The  time  response,  shown  in  Figures  11  and  12, 
is  markedly  improved  but  the  specifications  have  still  not 
been  met. 

For  the  eigenvalues  of  the  right-hand  column  of 
Table  VI,  the  state  weighting  matrices  have  been  increased 
by  10  again.  The  eigenvalues  are  all  stable  and  have  moved 
even  closer  to  the  origin.  They  are  again  plotted  in 
Figure  8.  The  time  response  for  this  system  is  presented 
in  Figures  13  and  14  and  shows  that  the  specifications 
have  been  approximately  met. 
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Table  VI 


Eigenvalues  of  Design  II 


Q  = 

1. 

0  I 

Q  = 

=  10.0  I 

Q  = 

qob 

»  1 

.0  I 

qob 

= 

10.0  I 

qob 

Eigenvalues 

of  $  +  r  f 

c  c 

.98892 

+ 

.11629i 

.98014 

+ 

. 11519i 

.95726 

.91629 

± 

. 33980i 

.87409 

+ 

. 31830i 

.76732 

.98613 

± 

. 14568i 

.97988 

+ 

. 14468i 

.96482 

.90177 

± 

.  36  4  3.0  i 

.85225 

+ 

.  33620i 

.73133 

Eigenvalues 

of  $  -  KC 

c  c 

.76513 

± 

.224861 

.30172 

+ 

Oi 

.04862 

.95803 

± 

. 11399i 

.93797 

+ 

.  12381i 

.93596 

.72822 

± 

. 21638i 

.86204 

+ 

Oi 

.88921 

.96553 

± 

. 14334i 

.95463 

+ 

.  14609i 

.95291 

.23944 

+ 

Oi 

.03449 

.85228 

+ 

Oi 

.87946 

Eigenvalues  of  4>g 

.95496 

± 

. 29170i 

.95496 

+ 

.  29170i 

.95496 

.84081 

± 

. 53607 i 

.84081 

+ 

.  53607 i 

.84081 

.86809 

± 

. 49121i 

.86809 

+ 

.  49121i 

.86809 

.83892 

± 

. 53899i 

.83892 

+ 

. 53899i 

.83892 

100.0  I 
100.0  I 

. 11396i 
. 22188i 
. 14330i 
. 21138i 


+  Oi 

±  . 12776i 
+  Oi 

±  . 14725i 
+  Oi 
+  Oi 


±  . 29170i 
±  . 53607i 
±  . 49121i 
±  . 53899i 


•  Q  =  qob  -  1 
o  Q  =  Q0B  =  10  1 

X  Q  =  Q0B  =  100  I 


Figure  8.  Z-Plane  Plot  of  Eigenvalues  of  Design  II 

48 


losx  *io' 
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VI  Design  III 


In  the  previous  two  designs  the  discrete  performance 
indices  have  been  of  the  form  of  the  simplified  performance 
index  of  Chapter  III,  Eq  (52).  This  approximation's  validity 
depends  on  the  accuracy  of  the  assumption  in  Eq  (47).  To 
see  the  effect  of  this  approximation,  Design  III  makes  use 
of  a  much  more  accurate  performance  index. 

The  form  of  the  discrete  performance  index  was  shown 
in  Chapter  III  to  be 

T={f  xT(kT)  xOcj)  +  xT0<~r)  Sj  M  Ckr) 

k=l  4  ~ 

4-  U&T)  s J  X  CWT)  +  Mr0<J)  R^aCkT)  0^) 

where 

sA  =  st  i  ct)  At 

Sj  =  Tl7?*)  a  ret)  At 

't> 

Rj  -  l^Ci)  a  P(P  +  R  Jt 

Equations  (143),  (144),  and  (145)  can  be  approximated  by 
numerical  integration  equations 

~  JZ  {lT(o)  Q  ±(o)  t  5  '  i~)  Q  Ktb')  +  .  .  ,  +  iT(^ )  Q  ICt?)]  ^0 


(Ml) 
(tv* ) 
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f 


Sj  —  ts  {iT(o)Q  rYo)  + 1  (ji) Q  V(-ii) +-  •  •  •  t-  r^)} 

Rj  « iS  frT< lo)  uRoif  iT(x)  $  r®  v  •  • .  +  rci?)  a rt£)}  +  kt 

This  performance  index  is  not  of  the  form  of  Eq  (52)  for 
which  the  solution  is  known.  However,  by  defining  a  matrix 
W  as 

=  Sj 

then,  to  complete  the  square  in  Eq  (142) 

Lw  xCkT)  t-  >u(kT)]  R  J  [  U  Ki\j)  +*  AX  (kT)]  -  /  (kT)  U  x  (kT) 

+-  a 0<T)  Rj  Wjf  CWTJ  -t-  x  CkT)  wTr<j  U 0<T)  t  xTllT)  \lj  U  (. kT) 

Now  the  performance  index  of  Eq  (142)  can  be  written 
I  ^CkT)^  -  W^U]  x  ( 1<T )  +■ 

CW  *0<T)  ^>uO<Tn'r^  tUxCUT)  +  4  CkT)] 

By  defining  a  new  control  variable 

f  CkT)  =  WxCkT)  +  A  CkT) 

the  performance  index  can  be  put  into  the  form  of  Eq  (52) 

Jz  i  f-  X^WT)  Q'xtkT)  +  ?TCk')Rj  €(kT)  ( 

k- 1 


0*0  ) 
(Wf) 

0*h) 

(iso) 

CisO 

(is  2) 

053) 

f  54  ) 
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whore 


Q'  -  cQa  -  (\A  U 

0  55 

This  performance  index  is  of  the  form  of  Eq  (52). 

The  solu- 

tion  for  a  state  equation  of  the  form  of  Eq  (37) 

is  given 

by 

H*(kT)  =  FxCkT) 

(  IT6) 

XJL*UT)  =  (F-w)xO<T) 

(157) 

where 

F  =  -  Ui+-  R/'r^pr1]  'r;1  rTp[t-  ru] 

USS) 

and  where  P  is  the  solution  to  the  disc  -etc  steady-state 
Ricatti  equation 

p  =.  [$- ru H - r wl - C« - p wlf r t  rr-nfVpEs-rulta'  osi) 

Design  III  begins  with  the  same  state  equation  as  the 
previous  two  designs 

XfcIa+OT]  =  fexcCkT)  f  Pc  -wCkT)  (UO' 

The  output  equation  is  the  same  as  for  Design  II 

yC  WT)  -  C,  x^Ofr)  t  Cs  XsCk T)  06  1 ) 

In  order  to  eliminate  the  control  spillover  the  same 
transformation  as  in  Design  II  is  made  to  the  control 
vector,  u 

AJ.  Ck T)  =  Alv(WT)  0;'-2N 
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I 


The  state  equation  becomes 

XjTO^0"r]-  ?c  Xc(WT)  -H  rc  M  v(kT)  063) 

and  the  performance  index  for  this  state  equation  is 

J  =  7  %  +  XcO<T)  v  CkV 

+-  vT^T)i  XXC6VT)  +  yT(kT)  R  V0<T)  OH) 

where 

(ibs) 

As  shown  earlier  in  this  chapter,  this  performance  index 
can  be  written 

J  -  *  H  /CkT)  Q'x  CkT)  +■  z'Oc'D  ^  (  14  ^ 

k  -  (  ~c 

where 

wr  =  S  R  068^ 

2(kT)-  U  X  ^T)  f  Y  (WT>  O^) 

and  the  state  equation  is  now 

Q.  l7  R^U  (no) 
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The  solution  is 


(177  ) 
(<12) 


?VOcT)  =  F  X  j:\c~T) 

V  *(W~0  -  CF-U)  X^kT) 

u*~Ck~r)  =  A(F-u)it.C  i<t) 

where 

F  =  -\j_+  ?v^n]  \^TrTr  [it-  rc^u]  o->?) 

and  where  P  is  the  solution  to  the  discrete  steady-state 
Ricatti  equation 

P  -  [i&-  rcMw]TP[.5c-  pca\u] 

-  \JC-  rc rn  (R^i-  ^iT rf ppm)’  '/F  r  r r  a i  >']  +  c{  m 6 

The  state  equation  can  now  be  written 

Xjtt.OT]  =  I<.XtGcT!  +  rrt(F-k)?,(lT) 

where  F  is  given  by  Eq  (173),  W  is  given  by  Eq  (168),  and 
M  is  given  by  Eq  (108)  from  Chapter  V.  The  observer  state 
vector  in  the  above  equation  is  given  by  the  observer  state 
equation 

iEiunr]  =-  Ic2ttkf)  C  utn n 

and  the  observer  output  equation  is 

y  (kT)  -  Cc  XC0<T) 

Substituting  Eqs  (161)  and  (179)  into  Eq  (178),  the  observer 
state  equation  can  be  written 


^ct/^OTj-=  C±c-KCc- QrtU  i- h)  x/),l)  +  kCcxc0<7)  f  KCS  x  0<t)  (iio)  j 


The  error  state  vector  is  again  defined  as 


e  am  =  xtam  -  xccut) 


(  hi) 


Then 


eZ(kii)T]~  xJ(kii)T]-  X<fr^«-0T1 

r  !cxcCk'T)f  r^rt(F-u)  £cCkT) 

"Tlc-  ^Cc  ^  /c(!<t)  -  KC6a«(I<t)  -)<£sxs(kT) 

-  Ot~  KCj  elVT)  -  Kc,  x,(kT)  (!t2) 

In  order  to  eliminate  the  observation  spillover,  the  gain 
matrix,  K,  is  again  defined  as  the  product  of  two  matrices 


K=  L/1 


where  MQB  is  as  given  in  Eq  (139)  in  Chapter  V.  Then  the 
error  state  equation  becomes 

eHtk.OT]=  Cfc-Kcj  etkT) 

As  in  the  previous  chapters,  this  equation  will  have  the 
same  closed-loop  eigenvalues  as  the  equation 

u]>4»t]=  Ij-  cJ^CkT) 


(if  3; 


ClH) 
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whore 


<j(kT)  -  (>‘6) 

With  K  defined  as  in  Eq  (183),  this  equation  can  be  trans¬ 
formed  into 

wlOc.o-r-]  -  l^uoro-  ucvr) 

where 


Un) 

Ott) 


The  performance  index  for  this  state  equation  is 
3  -  \  §  ur0J)  a  I  u>  (_I<T)  +-Hr0<~)  Srf  lnO<T) 

J;  |  ~  ~ 

i- blurts/  wturt'-t-  fa/urt  r„  Kart  mi) 

-  /'oft  —  "oft  - 

where 


Jl%t  "Jot 
R*.„  = 

As  shown  earlier,  this  perfoi’mance  index  can  be  written 

i  =  i  i  /art  <8  ytkT  1  *  §.>T)  ^  «rt 


(no) 

on) 


(HO 


where 

=  S„  Ou) 

''ob  '  <?& 

?  art  =  h,„  jsco<T)  +  uart  fu") 

®-rtS'Vw'-TV,Ht  on' 
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J 


b*(kT)  =  (kJ-  0.T)  (Mi) 

1  'aT 3  =  cIST-woa)y  aT)  on) 


where 

K  =  &  +  ^/i04ct  p«h;^„cjth  +  c^uJ  a-) 

and  where  PQB  is  the  solution  to  the  discrete  steady-state 
Ricatti  equation 

« 

-  =  [f:+  kT)] 


HI,  |<T+  <5' 

"o<j  ^oa 

The  matrix  L  of  Eq  (183)  is  thus  given  by 


0*0 


L=  (><"-  uo‘ 


'°a'  (j?o2) 

error  state  equation  is  now  as  given  in  Eq  (184) 
with  the  K  matrix  as  given  in  Eq  (183).  The  suppressed 
state  equation  is  now 

Xsfrk+0Tj=  5,K507)t  P^aT) 

~  ls  *s(.vr)  Ooi'j 
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i 


The  total  system  state  equation  can  now  be  written 


in  matrix  form  as 


^lOctOTl 

r*  t  n./K'F-U)  1  -rn(r^j!  © 

_ L 

— 

1  1 

O  |  1  ° 

i  cut) 

I  !  , 

L  °  i  0  i  \ 

_  yr). 

System  Performance 

In  Chapter  III  it  was  stated  that  the  simplified 
performance  index  was  a  good  approximation  so  iong-  as  the 
assumption  in  Eq  (47)  was  valid.  To  test  the  accuracy  of 
the  simplified  performance  index  of  Design  II,  Design  III 
was  run  at  the  same  sampling  rate  and  with  the  same  state 
weighting  matrices  as  those  of  Design  II  in  the  preceding 
chapter . 

The  largest  eigenvalue  of  Ac  is  -.01924  ±  3.84834i. 
With  a  sampling  time  of  0.1  seconds,  the  magnitude  of  the 
product  of  eigenvalue  and  sampling  time  is  equal  to  .38484. 

The  eigenvalues  of  Design  III  for  the  different 
weighting  matrices  are  shown  in  Table  VII  and  plotted  in 
Figure  15.  Comparison  with  the  corresponding  eigenvalues 
of  Design  II  in  Table  VI  shows  very  good  agreement  between 
the  two  designs  for  all  weighting  matrices.  Figures  16 
through  21  show  the  ;ime  response  of  Design  III  for  each 
different  Q  and  Qq^  .  In  each  case,  the  ccrresponding 
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Table  VII 


Eigenvalues  of  Design  III 


Q  =  1.0  I 

qob  -  1.0  I 


Q  =  10.0  I 

qob  =  io.o  i 


Q  =  100.0  I 
Qqb  =  100.0  I 


Eigenvalues  of  $  +  r  F 

c  c 


98892 

+ 

. 11629i 

.98014 

± 

.115191 

.95730 

± 

. 11401i 

91628 

± 

. 33983i 

.87392 

± 

. 31848i 

.76457 

± 

.  22166i 

98613 

± 

.  14568i 

.97988 

± 

. 14469i 

.96488 

± 

. 14333i 

90175 

± 

.364 34 i 

.85196 

± 

.  33647i 

.72697 

± 

.  21018i 

Eigenvalues 

of  $  -  KC 

c  c 

76125 

± 

. 23961i 

. 33956 

+ 

Oi 

.07775 

+ 

Oi 

96046 

± 

.114081 

.94487 

± 

.11869i 

.94238 

± 

.12067i 

72265 

± 

.239571 

.80540 

+ 

Oi 

.83848 

+ 

Oi 

96762 

± 

. 14348i 

.96003 

± 

. 14442i 

.95882 

.14480i 

. 28392 

+ 

Oi 

.06949 

+ 

Oi 

. 78030 

+ 

Oi 

.80943 

+ 

Oi 

Eigenvalues  of  <J>g 


.95496 

± 

. 29170i 

.95496 

+ 

.29170i 

.95496 

+ 

. 29170i 

.84081 

± 

. 53607i 

.84081 

4- 

. 53607 i 

.84081 

+ 

. 53607i 

.86809 

± 

. 49121i 

.86809 

± 

,49121i 

.86809 

+ 

.49121i 

.83892 

± 

. 53899i 

.83892 

+ 

. 53899i 

.83892 

+ 

. 53899i 
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Figure  16..  LOSX  VS.  TIKE,  Design  III,  Q  =  Q  =  I,  t  =  .1 
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Figure  19..  LOSY  VS.  TIME,  Design  III,  Q  =  Qqb  =  10  I,  T  =  .1 
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Figure  20.  I, CSX  VS.  TIME,  Design  III,  Q  =  Q  ^  =  100  I,  T  =  .1 
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VII  Effects  o f  Sampling  Time 


Aliasing 

The  size  of  the  sampling  time  plays  a  critical  role 
in  determining  system  response.  Theoretically,  the  lower 
limit  on  desired  sampling  time  is  T  =  0,  for  which  the  con¬ 
troller  becomes  a  continuous-time  system.  The  upper  limit 
on  desired  sampling  time  depends  on  the  natural  frequencies 
of  the  system  to  be  controlled.  To  illustrate  this  depen¬ 
dency  consider  a  continuous-time  system  with  a  time  response 
given  by 

X,0)  —  si*  (.*->,£  +■  ©)  (7°  ^ 

If  this  system  were  discretized  at  a  sampling  frequency, 

2t r  ws 

ojs  =  ,  such  that  0  -  wi  -  ~2  ’  't*le  discrete  time 

response  would  be 

X,  Cvro  =  *  '  ^  C  '<  T  +  G  )  C?o±) 

Next,  consider  a  continuous-time  system  given  by 

X2  it)  =  sin  HCw,  b  +  ©I  *  j  =  1, V,  ...  (20s') 

If  this  system  were  discretized  at  the  same  sampling  rate, 
the  discrete-time  response  would  be 

(kT)  -  s  im  -1-  l<~  +-  &  J 

—  Sim  [J  kT  u  ^  1-  ©  +  r  I r  i<  j  3 
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X  ( )<T )  -  s  .‘o  ClcT  u  +  &) 

7  1 


C  7°0 


Note  that  Eqs  (20G)  and  (208)  are  identical.  As  far  as  a 

digital  controller  with  the  same  sampling  rate  is  concerned, 

these  two  systems  are  the  same.  But  actually  they  are  quite 

different.  This  effect  is  known  as  aliasing.  For  an  open- 

loop  system  which  is  stable,  this  effect  will  degrade  the 

closed-loop  response  as  T  is  increased.  But  for  an  unstable 

open-loop  system,  this  effect  will  cause  the  closed-loop 

response  to  be  unstable.  Thus,  the  upper  limit  on  the 

desired  sampling  time  is  set  by  what  is  known  as  the  Shannon 

Sampling  Theorem.  Basically,  it  is  desired  to  always  choose 

a  sampling  frequency,  to  ,  such  that 

s 


where  <o  is  the  highest  natural  frequency  in  the  open-loop 
c 

system.  Therefore,  the  desired  sampling  time  must  be 

T  ^  TT  (7  t  o) 

I  —  £J 

c. 

In  the  following  analysis,  the  system  response  is 
examined  for  Designs  II  and  III  for  sampling  times  of  .1, 

.3,  .5,  and  .7  seconds.  In  all  cases,  the  state  weighting 

matrices  are  Q  =  Qqb  =  100  I.  For  this  model,  the  sampling 
time  of  .7  seconds  exceeds  the  condition  of  Eq  (210)  for 
modes  7,  8,  and  9. 
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System  Performance 


Presented  in  Table  VIII  are  the  eigenvalues  of 
Design  II  for  sampling  times  of  .3,  .5,  and  .7  seconds.  The 

plot  of  these  eigenvalues,  together  with  the  eigenvalues  in 
the  right-hand  column  of  Table  VI,  shows  the  general  trend  of 
eigenvalue  movement  as  T  increases.  Those  eigenvalues  that 
start  to  the  right  of  the  origin,  move  to  the  left  and  away 
from  the  real  axis  as  T  increases.  This  would  imply  that 
the  system  performance  is  being  degraded.  On  the  left  of  the 
origin  the  eigenvalues  continue  to  move  to  the  left  but 
begin  to  move  toward  the  real  axis  again.  For  eigenvalues 
on  this  side  of  the  origin,  this  also  implies  poorer  response 
However,  note  that  for  T  =  .7. seconds,  three  of  the  eigen¬ 
values  of  <b  are  to  the  left  of  the  origin  but  have  begun  to 
move  to  the  right  and  away  from  the  real  axis  indicating 
improving  system  response.  The  time  responses,  shown  in 
Figures  13,  14,  and  23  through  28,  reflect  these  results. 

As  shown  in  Chapter  V,  the  time  response  for  T  =  .1  meets 
the  desired  specifications.  For  T  =  .3  and  T  =  .5,  the  time 
response  is  progressively  worse.  However,  for  T  =  .7,  the 
time  response  is  progressively  worse.  However,  for  T  =  .7, 
the  time  response  is  actually  slightly  better  than  for 
T  =  .3. 

In  Table  IX  are  the  corresponding  eigenvalues  of 
Design  III.  Figure  29  plots  these  eigenvalues  and  shows 
the  same  trends  as  those  of  Design  II.  Figures  20  and  21 
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from  Chapter  VI  and  Figures  30  through  35  show  the  time 
responses  of  Design  III  for  the  above  sampling  times.  Again 
these  figures  show  the  same  trends  as  Design  II. 

Effects  of  Performance  Index  Approximation 

In  Chapter  III  it  was  shown  that  if  the  conditions 
of  Eq  (47)  are  met,  then  the  discrete  weighting  matrices 

,  and  ,  can  be  approximated  by  Eqs  (49)  through 
(51).  As  T  increases,  these  approximations  become  worse. 

To  examine  the  effect  of  the  increasing  disparity  between 
the  approximations  and  the  optimal  weighting  matrices,  com¬ 
pare  the  results  of  Design  II  and  Design  III.  Recall  from 
Chapter  VI  that  for  T  =  .1,  the  results  of  Design  II  and 
Design  III  were  almost  identical.  Now  compare  the  two  time 
responses  for  T  =  .3.  For  T  =  .3,  the  product  of  the  largest 
eigenvalue  with  the  sampling  time  is  equal  to  1.1545  which 
violates  the  condition  of  Eq  (47).  Thus,  the  approximations 
of  the  weighting  matrices  can  no  longer  be  considered  very 
accurate.  However,  in  examining  the  time  responses  it  can 
be  seen  that  they  are  still  quite  similar.  They  have  the 
same  general  form  and  differ  only  slightly  in  amplitude. 

Thus,  although  the  approximations  of  the  weighting  matrices 
have  declined  in  accuracy,  they  are  still  adequate  to  pro¬ 
vide  system  performance  approximately  equal  to  optimal  per¬ 
formance  index. 
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Table  VIII 


Eigenvalues  of  Design  II 


T 

= 

.3 

T  = 

=  .5 

r 

r  = 

=  .7 

Eigenvalues 

of  +  r  f 

c  c 

83954 

+ 

. 31080i 

.69058  ± 

.  46021i 

.52615 

± 

.  56060i 

35627 

+ 

.  41297i 

.04233  ± 

.  43947i 

.25409 

± 

. 37689i 

31332 

+ 

. 37699i 

.00962  ± 

.40257 i 

.43217 

± 

.  70772i 

83790 

+ 

.  39627 i 

.65179  ± 

.  58802i 

.32576 

± 

. 32880i 

Eigenvalues 

of  $  -  KC 

c  c 

11587 

+ 

Oi 

.71273  ± 

.  47869i 

.58061 

± 

. 62675i 

81828 

± 

. 31328i 

-.09740  ± 

.  18358 i 

.05999 

+ 

Oi 

37675 

+ 

Oi 

.67699  ± 

.612381 

.71708 

+ 

Oi 

83436 

± 

.  39749i 

-.16727  ± 

.  04321i 

.47035 

+ 

.  78033i 

10188 

+ 

Oi 

- 

.03550 

± 

Oi 

29535 

+ 

Oi 

- 

.84339 

+ 

Oi 

Eigenvalues  of  3> 

s 

62710 

+ 

.  77323i 

.08774  ± 

.  98873i 

.47830 

+ 

.  86642i 

13046 

± 

. 98290i 

-.94080  ± 

.  29482i 

.66052 

.  72440i 

02579 

± 

. 99197i 

-.83277  ± 

. 53016i 

.87876 

+ 

.  43862i 

14073 

± 

. 98143i 

-.94572  ± 

.  27833i 

.64255 

+ 

.  74023i 

77 


1.0. 9  .8  .7  .6  .5  .4  .3  .2  .1  0.0. 1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 


•  T  =  .1 
»T  =  .3 
XT  =  .5 
®T  =  .7 


I 


Figure  22.  Eigenvalue  Movement  as  T  Increases,  Design  II 
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Figure  23.  LOSX  VS.  TIME,  Design  II,  Q  =  Q()  =  100  I,  T  =  .3 
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Table  IX 


Eigenvalues  of  Design  III 


T 

= 

.3 

T  = 

=  .5 

1 

t  _ 

=  .7 

Eigenvalues 

of  4  +  r  F 

c  c 

83972 

+ 

. 31154i 

.68960  ± 

.  46089i 

.52487 

+ 

.564631 

33978 

± 

.  40202i 

.04274  ± 

.  43852i 

.  25801 

+ 

.  31201i 

29519 

± 

.  36135i 

.00963  ± 

.  40266i 

.43282 

+ 

.  71685i 

83826 

± 

.  39684i 

.65120  ± 

. 58845i 

.31988 

± 

.  06738i 

Eigenvalues 

of  $  -  KC 

c  c 

07285 

± 

.  45029i 

.76801  ± 

.  51209i 

.62671 

+ 

. 67532i 

86933 

± 

. 32156 i 

-.42047  ± 

.  .  44190i 

.52613 

+ 

Oi 

00061 

± 

.  47548i 

.70747  ± 

.  63890i 

.71005 

+ 

Oi 

86448 

± 

.  40809i 

-.51000  ± 

.  36803i 

.49143 

+ 

.  81507 i 

- 

.46452 

+ 

Oi 

- 

.78964 

+ 

Oi 

Eigenvalues  of  4g 

62710 

± 

.  77323i 

.08774  ± 

. 98873i 

.47830 

± 

.  86642i 

13046 

± 

.  98290i 

-.94080  ± 

.  29482i 

.66052 

± 

.  72440i 

02579 

± 

.  99197 i 

-.83277  ± 

.  53016i 

.87876 

± 

,43862i 

14073 

± 

.  98143i 

-.94572  ± 

.  27 833i 

.64255 

± 

.  74023i 
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Figure  34.  LOSX  VS.  TIME,  Design  III,  Q  =  0  =  100  I,  T  -  .7 
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Figure  35.,  LOSY  VS.  TIRE,  Design  III,  Q  =  Q0B  =  100  I,  T  =  .7 
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Comparing  the  time  responses  of  T  =  .5  shows  that 
for  this  sampling  time  some  differences  begin  to  appear  in 
the  two  design  performances.  For  T  =  .5,  the  eigenvalue¬ 
sampling  time  product  is  1.9242.  Finally,  for  T  =  .7,  the 
two  time  responses  show  significant  differences.  The 
eigenvalue-sampling  time  product  for  this  case  is  2.6939. 
From  these  results  it  would  seem  that,  though  the  weighting 
matrix  approximations  become  poor  as  AmaxT  approaches  1, 
they  are  still  good  enough,  even  slightly  above  1,  to  give 
system  performance  that  is  similar  to  that  of  a  design  using 
more  optimal  weighting  matrices. 


VIII  Conclusions 


Several  conclusions  can  be  drawn  from  the  preceding 
analysis.  First,  as  was  the  case  in  continuous-time,  the 
effects  of  control  and  observation  spillover  are  destabi¬ 
lizing.  For  a  fixed  number  of  sensors  and  actuators  with 
fixed  orientation,  the  singular  value  decomposition  trans¬ 
formation  can  completely  eliminate  these  effects  provided 
the  system  is  fully  controllable.  Elimination  of  the  spill¬ 
over  terms  guarantees  closed-loop  stability  of  a  system  which 
has  open-loop  stability. 

For  an  unstable  open-loop  system,  the  spillover  terms 
are  not  the  only  destabilizing  factors.  A  large  sampling 
time  may  drive  the  closed-loop  system  unstable  even  if  the 
spillover  terms  are  eliminated.  Even  for  stable  open-loop 
systems,  a  large  sampling  time  will  have  a  degrading  effect 
on  closed-loop  performance.  To  avoid  this,  a  sampling  time 
must  be  chosen  which,  at  the  very  least,  obeys  the  Shannon 
Sampling  Theorem.  Preferably,  a  sampling  time  well  within 
this  criterion  should  be  used  to  give  the  best  system 
performance  possible. 

Finally,  the  simplified  discrete  performance  index 
proved  to  be  very  accurate  for  the  case  where  the  sample 
time  is  small  compared  to  the  characteristic  times  of  the 
system.  Even  for  cases  where  this  condition  was  slightly 
exceeded,  the  approximation  yielded  performance  very  close 
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1" 


to  optimal.  In  general,  the  performance  was  accurate  for 

Am„„T  51  1  •  Therefore,  if  the  product  of  the  largest 
max 

eigenvalue  and  the  sampling  time  is  less  than  one,  the 
simplified  performance  index  may  be  used  without  appreciable 
loss  of  performance. 
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APPENDIX  A 

Eigenvector  Results  of  NASTRAN  Analysis 
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Eigenvalue  1 


Eigenvector  1 


Eigenvalue  3 


Eigenvector  3 


Eigenvalue  5 


Eigenvector  5 


Eigenvalues  and  Eigenvectors 


1.370E+00  Eigenvalue  2  =  2.151E+00 


-2.4?1E-01 

3.999E-01 

4.279E-02 

2.309E-01 

1.452E-06 

-1.489E-01 

-1.963E-02 

8.329E-02 

3. 398E-02 

4.808E-02 

-7.213S-02 

Eigenvector  2  = 

6. 8132-02 

-3. 696E-02 

7. 000E-02 

4.347E-02 

2 . 253E-02 

4.397E-0 2 

-4. 721E-02 

-1.962E-02 

5.451E-02 

5.296E-02 

4.936E-02 

4.397E-02 

-4.722E-02 

8.789E+00 

Eigenvalue  4=1 

.266E+00 

'  6.368E-02' 

~2.746e-02 

3.678E-02 

-4.7582-02 

4, 000E-01 

-2.2492-05 

1.984E-01 

-1.718E-01 

1.145E-01 

2.977E-01 

2. 010E-01 

-6.8172-05 

1.548E-01 

-2.612E-C1 

6, 804e~02 

3.4362-01 

9.782E-02 

-8.190E-02 

1.3632-01 

-1.7162-01 

1 , OOOE-Ol 

3.S94E-01 

9.784E-02 

8.192E-02 

1 , 481E+01 

Eigenvalue  6=2 

.652E+01 

-8.783E-02" 

‘1.353E-05 

-5.070S-02 

1.218E-11 

-1.299E-01 

3.4022-11 

3.095E-01 

-2.0412-01 

1 . 786E-01 

3.535E-01 

-  -3. 5142-01 

-6.0572-06 

2.866E-01 

-2.041E-01 

1 .224E-01 

-3.535S-C1 

1.139E-02 

1.086E-04 

2.494E-01 

4.0822-01 

1.868E-01 

6.802E-10 

1.140E-02 

5.065E-10 
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Eigenvalue  7  =  3.222E+01 

-2.661E-02 
4.607E-02 
3.302E-05 
3. 374E-02 
-5. 844E-02 
Eigenvector  7  -  3.2312-05 

2.733E-02 

-5.481E-02 

-4.913E-01 

3.382E-02 

-5.108E-02 

4.909E-01 


Eigenvalue  9  =  7.992E+01 

~9.907E-02 
5. 720E-02 
1.729E-01 
1.076E-01 
6.213E-02 
Eigenvector  9  =  -4.953E-01 

-1.679E-01 
-2.1982-01 
-1 . 110E-02 
-2.743E-01 
-3.55^-2-02 
-1 ,1092-02 


Eigenvalue  11  =  1,1932+02 

'6.370E-02' 
3.678E-02 
9.588E-02 
-2.401E-01 
-1. 3862-01 
Eigenvector  11  =  -2 . 6052-01 

-8.6062-02 

3.9442-01 

6.970E-03 

2.9842-01 

-2.7192-01 

6.971E-03 


Eigenvalue  8  =  3.2612+01 

-2.994E-02 
-1.731E-02 
8.784E-02 
4.071E-02 
2.360E-02 
Eigenvector  8  =  3.554E-02 

2.7422-02 
2.7982-02 
-4.875E-01 
3.799E-02 
9. 8102-03 
-4. 879E-01 


Eigenvalue  10  =  1.062E+02 

-3.3902-03 

5.850E-03 

-1.6052-05 

-2.286E-01 

3.960E-01 

Eigenvector  10  =  4.9642-05 

3.7832-01 

4.5542-02 

-1.471E-02 

-2.2862-01 

-3.0492-01 

1.4722-02 


Eigenvalue  12  =  1.9512+02 

~3.206E-0f 
1.851S-02 
6.438E-02 
-4.0262-01 
-2,3242-01 
-1.3052-01 
3.204E-01 
-1.587E-01 
-9.2782-03 
2.272E-02 
3. 5682-01 
-9.282E-03 
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APPENDIX  B 

Matrices  B  ,  B  ,  B 
c  s’  r 
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\  THE  8  SUPPRESS E 8  MfcTklX  T5 


I 


Discrete  Cnen-Ioop  Eigenvalues 


Eigenvalues  of  S’ 

.99257  ♦  .ll6?3i 
.93571  +  .3^7691 
.98854  +  . l4604i 
.92508  +  .374681 

.93732  +  .343411 
.47991  +  . 871 24i 
.90275  1  .425031 
.40205  t  -909331 

.‘83110  +  .550851 
-.20473  +  .969771 
.74020  +  .666941 
-.34275  +  .929231 

.67984  +  .727781 
-.78551  +  ,59864i 
.51^97  +  .851231. 
-.88936  +  .427151 


Eigenvalues  of  3f_ 

_ _ s 

.1 

.95496  +  .291701 
.84081  +  .53607i 
.86809  +  .491211 
.83892  +  .538991 

.3 

.62710  +  .773231 
-.13046  +  .982901 
.02579  +  .991971 
-.14073  +  .981431 

.5 

.08774  +  .988731 
-.94080  +  .294821 
-.83277  +  .530161 
-.94572  +  .278331 

.7 

-.47830  +  .866421 
-.66052  +  ,72440i 
-.87876  +  .438621 
-.64255  +  .740231 
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